Superfluid phase diagrams of trapped Fermi gases with population imbalance 
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We present phase diagrams for population imbalanced, trapped Fermi superfiuids near unitarity. In addition 
to providing quantitative values for the superfluid transition temperature, the pairing onset temperature and 
the transition line (separating the Sarma and phase separation regimes), we study experimental signatures of 
these transitions based on density profiles and density differences at the center. Predictions on the BCS side of 
resonance show unexpected behavior, which should be searched for experimentally. 



PACS numbers: 03.75.Hh, 03.75.Ss, 74.20.-z 

The study of ultracold trapped Fermi gases, as they vary 
from BCS to Bose-Einstein condensation (BEC) is a rapidly 
exploding subject which is defining new directions in con- 
densed matter and atomic physics. With the recent focus on 
population imbalanced superfiuids JH 12, [H 01, this body of 
work has also captured the attention of physicists in other dis- 
ciplines [5]. Indeed, it is hard in recent years to find a subfield 
of physics which has such a broad appeal across the different 
sub-specialties. The goal of the present paper is to address the 
multiple superfluid and normal phases which appear in these 
trapped gases and which are viewed as possible prototypes for 
high temperature superconductors @], as well as quark and 
nuclear matter. The various phases we contemplate become 
stable or unstable as one alters the populations of the two 
different spin species or changes the temperature T. There 
has been a very extensive theoretical literature on this subject 
which is almost exclusively confined to T = 10, . Our 
emphasis here, is on finite temperature effects iot UOIl . which 
are essential in order to address the actual experimental situ- 
ation. We do this here along with presenting predictions for 
new experiments. 

The fermionic superfluid ground state wavefunction in the 
presence of BCS-BEC crossover effects (with population im- 
balance) is almost universally assumed 0, H 0] to be of the 
BCS-Leggett form. The excitations of this ground state are 
expected to be particularly interesting in the intermediate or 
unitary regime where they consist of non-condensed pairs as 
well as gapped fermions. While these finite momentum pair 
excitations are frequently omitted in the literature lfTll[T2tl . we 
stress that they play an essential role. Only with their inclu- 
sion will one find quantitatively meaningful estimates of T c ; 
this temperature is reduced dramatically (relative to a naive 
mean field estimate) as a result of the opening of an excitation 
or pairing gap well above the transition. Indeed, the separation 
of the pairing onset temperature T* and the phase coherence 
temperature T c is well documented in experiments which are 
consistent with a "pseudogap" (or normal state excitation gap) 
in the unitary trapped gases |@] . 

Without the contribution of non-condensed pairs J^,[l3t] one 
often finds bimodal particle distributions at p = 0. There- 
fore, theories which ignore these pairs fTTl Il2tl may not be 
applicable for establishing bimodal distributions associated 
with population imbalance. Pairing fluctuations were par- 
tially included (through the number equation) in Ref. fl4ll 
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based on the Nozieres-Schmitt-Rink (NSR) scheme lfl5ll but 
in the absence of a trap. However, the finite T NSR treat- 
ment is not designed to be consistent with the standard ground 
state If7l Isil [l6ll . In the presence of population imbalance, it is 
now well known that one has to consider homogeneous phases 
(which we refer to as the "Sarma") and Larkin-Ovchinnokov- 
Fulde-Ferrell (LOFF) phases [17] as well as phase separa- 
tion lfl8tl . In addition to these competing phases, the nor- 
mal phase may appear as a highly correlated or paired state 
without phase coherence, or as a simple Fermi gas (unpaired) 
phase. 

We begin with the (local) thermodynamic al potential (per 
unit volume) for the allowed states (excluding the more com- 
plicated LOFF phase which appears to be of less interest near 
unitarity and at all but the lowest temperatures [ 19])- We focus 
on a two-component Fermi gas in a harmonic trap potential. 
For a normal or superfluid phase in which pairing correlations 
are present the thermodynamic al potential (tltot) consists of 
two highly inter-dependent contributions: from the gapped 
fermions (£1 /) and the non-condensed pairs or bosons (f^) 

n tot = n f + n b (i) 
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Competing with this phase is the free Fermi 
phase which has thermodynamical potential £lf ree = 
— TJ2k a m (l + e~£ fc<T / T ). The (gapped) fermion and pah- 
dispersions are given respectively by E^ — \f%L + A 2 for a 
contact (short-range) pairing interaction with strength U and 
Cl q = q 2 /2M* — impair, where M* and fi pa i r are the effective 
mass and chemical potential of the pairs. Both M* and Z 
are obtainable from a microscopic T-matrix approach lllOll . 
Here E ka = =F h and £ kcr = £k T h for spin a =T, 4, 
respectively, where £k = k 2 /2m — fj,, fcg=?j=l with 

M = (MT + Mi)/ 2 and h = Ot - Ml)/2- 

We distinguish the order parameter A sc from the total gap 
A. Non-condensed pairs are associated with a pseudogap con- 
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tribution, A pg , to the total gap via A = A 2 
form of SI;, contains a free boson-like contribution; the pah- 
condensate does not contribute to fit, directly. Although it is 
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not a necessary assumption, in order to simplify the formal 
description we assume that M* depends on T only, as is rea- 
sonably consistent with our microscopic theory IU0I1 and fi pa ir 
is a function of A and T only. In the summation over boson 
momentum q, we will impose a cutoff q < q c , where q c is the 
minimum value of q which satisfies Eu +£k+q — — for at 
least one value of k; above this momentum, pairs become dif- 
fusive by decay into the particle-particle continuum. It should 
be noted that £1 f and fib couple only via fi pa ir- Above T c this 
leads to an extra contribution to the gap equation. 

Our self-consistent equations can be expressed as varia- 
tional conditions on QtoU which yield equations for the to- 
tal excitation gap A, the pseudogap A pg , the fermion number 
density n and population imbalance Sn, respectively: 

dfltot _ „ dfttot _ q dfltot _ _ dfttot _ ^ 
OA d/ipair <9/i dh 

Importantly, below T c we have /i pa ;r = 0, and above T c we 
have A = A ps . These equations reduce to the usual temper- 
ature dependent BCS-like equations for the gap and number 
equations in the literature. The effect of the pseudogap, which 
is new to the present formalism, has also been extensively dis- 
cussed in earlier papers |sll0|], although here we have pre- 
sented it in a different, but equivalent fashion. 

In a trap, T and h are independent of the radial position. 
The distribution of the local /i(r) is determined by the force 
balance equation, — Vp = n\7V ex t, where p — —Qtot is the 
pressure and V ex t = \mio 2 r 2 is the harmonic trap potential 
with mean angular frequency ui. Using the number equa- 
tion n = —dCltot/dfi, we obtain V/i(r) = —W ex t(r), or 
jtz(r) = - V ext (r), where p = £t(°) and Kxt(0) = 0. 
This shows that the force balance condition naturally leads to 
the usual local density approximation (LDA). To ensure mean- 
ingful comparisons between states under different conditions, 
we fix the total particle number N — J d 3 rn(r) and the num- 
ber difference SN. The Fermi energy Ep = Tp = k F /2m 
is defined as that of a non-interacting system with the same 
N = N~i + N L at the polarization p = SN/N = 0. We 
assume > N± so that h > 0. 

The physical state corresponds to the minimum of iltot and 

free- At a particular trap radius when (and if) these become 
equal, the system will have a first order transition from a state 
in which there is pairing (but not necessarily superfluidity) to 
an unpaired Fermi gas phase. We assume an infinitely thin 
interface (i.e., domain wall) for this phase separation. There 
is as yet no complete microscopic theory for the interface en- 
ergy, so that we do not include it here. The philosophy behind 
our treatment of phase separation is very similar to that in pre- 
vious T = papers Q8L |20J] except that we include the impor- 
tant effects of temperature (in a fashion consistent with the 
well studied BCS-Leggett ground state). Across the interface, 
we require thermal, chemical, and mechanical equilibrium so 
that T, [i a and p are continuous. 

Fig.Q]shows the phase diagram at unitarity for a population- 
imbalanced Fermi gas in a trap. Here 1/kpa = 0, where a is 
the s-wave scattering length between fermions. Phase sepa- 
ration (labeled PS) occupies the lower T portion of the phase 
diagram, where the gap A jumps abruptly to zero at some trap 
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Figure 1 : (Color online) Phase diagram of a population-imbalanced 
Fermi gas in a harmonic trap at unitarity. The solid lines separate dif- 
ferent phases. Above the (red) dashed line but within the PS phase, 
the superfluid core does not touch the domain wall. Here "PG" indi- 
cates the pseudogapped normal phase. The black dot labeled "TCP" 
indicates the tricritical point. The four points indicated by the trian- 
gles labeled (al)-(b2) correspond to the density profiles in Fig. [2] 



radius. At intermediate T, there is a (yellow-shaded) Sarma 
phase, where A vanishes continuously within the trap. It 
evolves into a (dotted) pseudogap (PG) phase as the superfluid 
core vanishes at higher T. A normal (N) phase without pairing 
always exists at even higher T. In contrast with earlier work 
OH Hi, here we include pair fluctuations or non-condensed 
pairs, which are essential in order to obtain quantitatively cor- 
rect values for T c as well as for the finite temperature particle 
density profiles Ball 311. 

We stress that with or without the trap, the PS phase is the 
ground state at unitarity for any polarization. The differences 
between the trap and homogeneous phase diagrams are prin- 
cipally quantitative. The finite temperature boundary between 
the Sarma and PS phase is shifted as a result of the trap poten- 
tial, allowing the Sarma state to exist at low T and low polar- 
izations. The location of the tricritical point (TCP) also shifts 
from low polarization in the uniform case (not shown) to high 
polarization in the trapped case. Away from p = there is a 
minimum temperature required to arrive at this Sarma phase, 
as a consequence of the well documented negative T = su- 
perfluid density. 

In Fig. |2] we present four representative density profiles 
which show the behavior on different sides of the transition 
line between the PS and Sarma phases. The (al)-(a2) pah- 
have the same T but different p, while the (bl)-(b2) pair have 
the same p but different T. We have chosen these points some 
distance from the transition line in order to illustrate the dif- 
ferences. It should be noted, however, as points (al) and (bl) 
(which are within the Sarma phase), move closer to the tran- 
sition line, A drops rapidly (but not discontinuously), as a 
precursor to abrupt phase separation. The transition from the 
pure Sarma phase to the PS phase is, thus, a relatively smooth 
one. These LDA-based calculations should apply to situations 
where the trap geometry is reasonably isotropic, as for exam- 
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Figure 2: (Color online) Three-dimensional density profiles corre- 
sponding to the four points labeled in the phase diagram at unitarity 
shown in Fig. [TJ (al) T = 0.1T F , p = 0.5. (a2) T = 0.1T F , 
p = 0.7. (bl) T = 0.12Tf, p = 0.6. (b2) T = 0.08T F , p = 0.6. 
The (black) solid, (red) dashed, and (blue) dot-dashed lines corre- 
spond to rif, nj,, and <5n, respectively. Here Rtf = \/1Ef /muj 2 
is the Thomas-Fermi radius, and the units for density are kp 3 . 



pie in the MIT experiments dUt] and where one might be able 
to ignore surface energy contributions. On the other hand, this 
behavior appears at odds with recent experiments from Rice 
[2J, |4J which report pronounced changes in the aspect ratio 
of the profile as the transition line is crossed. These changes 
have been attributed to extreme trap anisotropy and associated 
interface energy effects 12 111 . 

In a different class of experiments it was proposed that by 
measuring densities at the center of the trap as one sweeps 
T or p, one can infer when the Sarma-PS transition line is 
crossed [ 22] as well as where superfluid transition [i.e., T c (p)] 
occurs [4]. As in these experiments, we plot the ratio rif / at 
the trap center as a function of T and p in Fig. [3] As shown in 
Fig. [3^, when p is fixed at a relatively high value while sweep- 
ing T, the curve starts at 1 in the PS regime at low T, and be- 
gins to increase when the the transition line to the Sarma phase 
is crossed. This behavior reflects the fact that Sarma states can 
accommodate higher core polarizations than their PS counter- 
parts. In Fig. [3J3, when instead T is fixed below TCP while 
sweeping p, the curve is a straight line independent of p, as 
observed experimentally |4J]. Here one can infer, that the core 
has equal population in both the Sarma state at low T and the 
PS phase at any (T,p). This figure also underlines the fact 
that PS states exist up to very high polarizations (p > 0.99). 
Finally, when T/T F is higher than that of the TCP (Fig. 0;), 
the ratio increases rapidly after the superfluidity disappears, 
as appears to be observed experimentally 0]. Figure [5}; re- 
veals that when superfluidity is present, it resides in the center 
of the trap, expelling the excess fermions outside the core. It 
should be noted that if one performs the same experiment as 
in Fig. [3^, but with much lower polarization (not shown), we 
find that the crossing of the transition line will be barely ob- 
servable, as is necessary in order to be consistent with Fig.[3j3. 
No matter how the Sarma-PS transition line is crossed, at low 
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Figure 3: (Color online) Behavior of /nj at the trap center at uni- 
tarity, as a function (a) of T/Tf at fixed p — 0.5, of p at fixed (b) 
T = 0.057V and (c) T = 0.2T F . The (blue) dashed line in (c) 
indicates n-\/ny = 1. The arrows labeled "Transition line" indi- 
cate where the PS-Sarma transition occurs, and the arrows labeled 
"A sc = 0" indicate where superfluid condensate disappears in the 
Sarma state. 



p (and thus low T), the ratio n^jn^ remains 1. Therefore, this 
class of experiments may not map out the transition curve for 
all p. 

Our theory can be generalized to address the whole of BEC- 
BCS crossover. As one passes from unitarity towards the 
BEC regime we find the fraction of the PS phase in the phase 
diagram decreases progressively, disappearing at 1/kpa ~ 
2.04. For stronger couplings, the superfluid contribution to 
the phase diagram consists only of the Sarma phase, at all T. 
The same observation at T — was first reported in Ref. 12311 . 
This can also be inferred from the homogeneous phase dia- 
gram in Ref. j24ll . which shows that the Sarma state is stable 
at low T for any p providing 1/kpa > 2.3. 

In Fig. |4] we present the phase diagram for a population- 
imbalanced trapped gas on the BCS side of resonance, where 
l/fepo = —0.5. It differs significantly from the unitary case 
and is, in many ways even more rich. Importantly, at T = 0, 
the PS phase is no longer stable at very high p. We understand 
this by noting that in the BCS regime, the pairing interaction 
is relatively weak and the gap A small. At sufficiently high p, 
we have h > A(r = 0) so that an unpolarized BCS superfluid 
core can no longer be sustained. 

Equally important is the fact that in this BCS case, the 
Sarma and PS phases do not connect and an intermediate 
phase appears between the Sarma and PS states. As a conse- 
quence, the boundaries of the Sarma, PS, and (pseudogapped 
or unpaired) normal phases do not meet except possibly at 
p = T = 0. We presume that this intermediate phase is a 
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Figure 4: (Color online) Phase diagram of a polarized Fermi gas, as 
in Fig.[TJexcept that l/kF<i = —0.5. Here the Sarma and PS phases 
are separated by an intermediate normal regime. The (green) dotted 
line indicates a sweep of T at p — 0.5, with five possible structures 
labeled: (1) PS, (2) N, (3) Sarma, (4) PG, (5) N. Shown in the inset 
is n j /nj at the trap center as a function of T/Tf- 

normal Fermi gas (N). To understand its appearance, we note 
that as the BCS regime is approached, (i) phase separation be- 
comes problematic because finite temperature (which enters 
via T I A) has a stronger effect, allowing the polarization to 
penetrate into the center of the core and thereby making PS 
more difficult, (ii) In addition, the intermediate temperature 
Sarma phase becomes more fragile as the pairing is weakened. 



As a result the Sarma and PS states retreat from each other as 
seen in Fig. [4] We cannot rule out LOFF-like states as an al- 
ternative candidate in place of N. However, we have found 
B19I1 these phases in general have very low T c and should not 
persist at these higher temperatures. In the inset of Fig. [4] we 
plot n^/ni at the trap center as a function of T/Tp. The N 
state between the PS and Sarma phases would be manifested 
by sudden jumps at the PS-N and N-Sarma boundaries. This 
prediction can be used to test the existence of the intermediate 
N state in the phase diagram. 

We end with another prediction concerning how the "mixed 
normal" region of the trap, emphasized experimentally in Ref. 
yfl, evolves with T. As noted in earlier work [25], within the 
Sarma phase, this mixed normal state consists of highly corre- 
lated non-condensed pairs which necessarily become less sig- 
nificant as T decreases. However, with decreasing T, as seen 
from Figs.[T]and|4j the Sarma phase gives way to stable phase 
separation. Only in the sliver above the (red) dashed line will a 
phase separated gas contain a correlated mixed normal region 
(appearing between the superfluid core and the domain wall). 
Everywhere in the PS phase, there is a mixed normal region in 
the gas outside the domain wall, with no pairing correlations, 
as appears at strictly T = 0|1§]. Thus a change should oc- 
cur from a highly correlated to an uncorrelated mixed normal 
region at large radii as T is progressively decreased. Further 
experiments should help to address these predictions. 
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